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ABSTRACT 

We report on high resolution CO(l-O), CS(2-1) and Smrn continuum Combined Array for Research in 
Millimeter Astronomy (CARMA) observations of the molecular outflow host and nearest quasar Markarian 
231. We use the CS(2-1) measurements to derive a dense gas mass within Mrk231 of 1.8 ± 0.3 x 10^° Mq, 
consistent with previous measurements. The CS(2-1) data also seem to indicate that the molecular disk of 
Mrk231 is forming stars at about normal efficiency. The high resolution CARMA observations were able 
to resolve the CO(l-O) outflow into two distinct lobes, allowing for a size estimate to be made and further 
constraining the molecular outflow dynamical time, further constraining the molecular gas escape rate. We And 
that 15% of the molecular gas within the Mrk231 outflow actually exceeds the escape velocity in the central 
kiloparsec. Assuming that molecular gas is not constantly being accelerated, we And the depletion timescale 
of molecular gas in Mrk 231 to be 49 Myr, rather than 32 Myr, more consistent with the poststarburst stellar 
population observed in the system. 

Subject headings: galaxies: active - galaxies; evolution-galaxies; ISM - galaxies; nuclei - galaxies; quasars; 
individual (Markarian 231) 


1. INTRODUCTION 

Understanding the role of active galactic nuclei (AGN) in 
its host galaxy’s transition from a blue starforming disk to 
a red quiescent ellipsoid is one of the forefront problems in 
studying galaxy evolution, with molecular gas at the center 
of the debate. A dearth of galaxies appearing in the opti¬ 
cal “green valley” (lEaber et al.ll2007^ seems to indicate that 
galaxies must transition between these two stages rapidly. It 
has been suggested that the expulsion of a galaxy’s starform¬ 
ing material through AGN feedback might provide a mech¬ 
anism that is both sufficiently powerful and rapid to ac hieve 
this transition (iSpringel et al.l20^lHopkins et al.ll2006h . Re¬ 
cent observations ha ve shown many galaxies with AGNs host 
molecular outflows (lEerug lio et al l2010t lAlatalo et al]l201 It 


lAalto et ^l2012bl;ICicone et al.ll2014l and references therein). 

In fact, AGN-driven molecular winds might be a ubiquitous 
property, as molecular observations of these types of objects 

become more numerous. _ 

Ma rkarian 231 (at a distance of 181 Mpc; iBoksenberg et al.l 
is an ideal laboratory to study the ways in which 
supermassive black holes and star formation interact with 
molecular gas. Mrk231 has bo th a rapidly accreting, 
powerful AGN (iBraito et al.l 120041) and a starburst, with 
a sta r formation rate of 170 Mq yr“^ (iVeilleux et al.l 
IM)1 . and is classified an ultral uminous infrared galaxy 
(ULI RG; ISanders & Mirab ell 19961) . Its disturbed morphol¬ 
ogy (iHutchings & Nefll 11987 ) conf irm that it is in the late 
stages of an ongoing merger, which [Hamilton & Keell (Il987h 
suggest is the origin of both the starburst and quasar activ¬ 
ity. Recently, it has also been shown that Mrk 231 hosts 
a massive outflow, observed in molecular gas (lEischer et al.l 
l2010t lEeruglio et^ l2010t lAalto et al.l l2012al). neutral gas 


(iRupke et al.l2005l;lRupke & Veille uxl201 ll Teng et al.l2oT3l) 

and warm ionized gas (iLipari et al.l 120091 iRupke & VeilleuxI 
120111; IVeilleux et al.ll2013l) . These properties make Mrk 231 
the ideal laboratory to test the ways in which dense molecular 
gas, star formation, and an AGN intermingle within a merging 


system, shedding light on the fate of the starforming material, 
and ultimately, the galaxy as a whole, after a merger. 

We present new high resolution Combi ned Array for Re- 
search in Millimeter Astronomy (CARMA; iBock et al.ll2006l) 
continuum, CO(l-O) and CS(2-1) observations of Mrk231. 
In |j2] we describe the observations and data analysis. In 
we use the CARMA measurements to put the molecular gas 
in Mrk 231 into context, both in terms of star formation prop¬ 
erties and depletion timescales. In SjH we present our conclu¬ 
sions. We use the cosmo logical parameters H o = 70 km s“^, 
Om = 0.3 and A = 0.7 (ISpergel et al.ll2007l) throughout. 

2. OBSERVATIONS AND DATA REDUCTION 

Mrk 231 was observed with CARMA between 2010 
November and 2011 May in two different configurations; C 
array (1" resolution) and B array (O'.'?). ^^CO(l-O) and CS(2- 
1) were observed simultaneously, utilizing the upgraded cor¬ 
relator. The primary beam has a diameter of 2' at CO(l-O), 
which covers all emission from Mrk 231. All observations 
used a long integration on a bright quasar to calibrate the 
passband, and alternated integrations between a gain calibra¬ 
tor (1419-1-543) and Mrk231. We then used the Mrk231 3mm 
continuum source (iJovce et al.lll975[) to self calibrate. The 
self calibration step was able to mitigate atmospheric fluctu¬ 
ations for the longer baseline arrays. The data were reduced 
using the Multichannel Image Reconstructi on Image Analysi s 
and Display (MIRIAD) software package (iSault et al.lll995h . 
C alibration and data re duction steps were followed identically 
to lAlatalo et al.l (120131) . Eigure[T] presents the CO(l-O) spec¬ 
trum, using the aperture defined by the integrated intensity 
(momentO) map (seen as an inset in Eig. [T]). Table[T]presents 
the derived properties from the CARMA observations. 

The MIRIAD task uvlin was used to separate out contin¬ 
uum emission from Mrk 231 from the line emission, and we 
estimate the continuum flux of Mrk 231 to be 20.9 ± 0.3 mJjOl 
in the unself-calibrated data. A centroid was computed to de- 

* This does not include the absolute flux calibration uncertainty of 20% 
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TABLE 1 

Mrk231 properties 


Tracer 

^maj X ^min 
(") 

K per Jy^ 

Av 

(km s“^) 

RMS 

(mJy beam~^) 

RMS per channel 
(mJy) 

Area 

n" 

-^peak 

(Jy bearn”^ km s“^) 

-^int 

(Jy km s“^) 

CO(l-O) 

1.80 X 1.57 

35.54 

100 

2.77 

7.70 

1.1 

30.4 ± 5.8 

74.3 ± 3.0 

CS(2-1) 

1.40 X 1.26 

78.35 

34 

1.28 

1.82 

3.9 

0.92 ±0.13 

1.43 ±0.19 

Blue 

1.19 X 1.07 

78.48 

400 

0.57 

- 

5.3 

1.91 ±0.23 

4.30 ±0.53 

Red 

1.19 X 1.07 

78.48 

400 

0.54 

- 

5.9 

1.50 ±0.21 

4.60 ±0.53 


i Kelvin per Jansky factor for each observation 


termine the position of the 3mm c ontinuum and c ompare it to 
the VLBI-determined radio core (iMa et al.lll99i^ . and found 
that the two positions agree to 0.04". 

To image the blue-shifted and red-shifted wings, we created 
images with velocities between —800 < v < —400km s“^ 
and 400 < v < 800km s“^, respectively. Figure [1] shows 
the continuum-subtracted CO( 1-0) spectrum derived from the 
100km s“^ CARMA channel maps, using the momentO map 
(Fig. [T] inset) as the masking aperture in each channel. The 
corresponding wing images and spectra are shown i n FigurelH 
overl aid on B-band HST imaging of Mrk 231 (right: iKm et al.l 
l2013h . CARMA derives a total flux of 74.3 ± 3.0 Jy km s“^, 
in agree ment within errors with the flux recovere d by Plateau 
de Bure (iFeruglio et akllMTlil: ICicone et al.ll2012h . 

CARMA has successfully resolved the high velocity CO(l- 
0) wings, finding a separation between the centroids of 0.49" 
(415 pc), a distance 5.8 times larger than the centroiding er- 
ro]0, and the blue- and red-shifted lobes have sizes of 2.3" 
(1980pc) and 2.2" (1820pc), respectively. The fluxes asso¬ 
ciated with each imaged wing was determined by summing 
the emission that lay inside the 3a contours of the emission, 
and were found to be Fqo .blue = 4.30 ± 0.53 Jy km s ^and 
FcO; red = 4. 60 ± 0.53 Jy km s which is larger than what 
was found in ICicone et al.l (120121) . though it is possible that 
our apertures were larger. The lobes also appear to be sep¬ 
arated in the East-West direction, with the centroid of the 
blue lobe located at [12:56:14.245, -1-56:52:25.24] and for the 
red lobe centroid at [12:56:14.213, -1-56:52:25.21], with the 
3mm continuum point source (white cross) located between 
the lobes. 

CS(2-1) observations are shown in Figure |3] Construc¬ 
tion of channel maps, integrated intensity (momentO) maps, 
mean velocity (momentl) maps, spectra, and the correspond¬ 
ing determinations of the root mean square no ise were per¬ 
forme d identically to the methods described in lAlatalo et al.l 
(120131) . CS(2-1) channels have an (optically-defined) veloc¬ 
ity width of 34km s“^. The total CS(2-1) flux, determined 
by summing across the shaded channels in the spectrum, is 
1.43 ± 0.19 Jy km s“^, and covers an area of 3.9 arcsec^. 

3. RESULTS AND DISCUSSION 
3.1. Dense molecular gas traced by CS(2—1) 

Using the local t hermal equilibrium f rsi'b-u-ATHo.riRn.ss 
relation discussed in lAlatalo et al.l (120151) of 

Ncs = 1.90 X 10“Tex /cs( 2 -i) 

where /cs( 2 -i) is the integrated intensity of the CS(2-1) line 
in K km s^“t^ assuming Tex ~ 70K (Ivan der Werf et al] 
120101) . and a CS abundance of Ri 10“®, we And a total dense 
gas mass of ^ 1.8 ± 0.3 x 10^° Mq. This dense gas mass 

^centroid “ ^beam/(2 X SNR.) 


is consi stent with the total dense gas mass found with other 
tracers (iS olornon et al.lll99^ iSolomon & Vanden Boutil2005l: 

IFeruglio et al.l 120101) . given the uncertainty in the CS/H 2 
abundance, which is at least a factor of 2. Our dense gas 
mass estimate is lower than the dynamical mass derived 
within 1100pc of the cente r of Mrk231 is 3.15 x 10^° MgO 
(iDownes & Solomonll998h . necessitating that the majority of 
the central mass in Mrk 231 is in the form of dense molecular 
gas. The signal-to-noise ratio of the CS data and resolution 
of these data are not able to confirm whether the CS is p art of 
the face-on wa rped disk reported i n iDavies et al.l (120041) . and 
in dense gas by I Aalto et al.l (120151) . Despite the limitation in 
resolution, the gradient seen in the CS momentl map also ap¬ 
pears co nsistent with the HCN kinematics seen bv lAalto et al.l 
(I2012al) . The agreement of the molecular gas mass also ap¬ 
pears consistent with the idea that the majority of the molec¬ 
ular gas in Mrk 231 is in a dense form. 

Our integrations were not sufficiently sensitive to detect 
CS(2-1) emission associated with broad wings (as is the 
case in other dens e gas tracers, including HCN, reported by 
I Aalto et ^I2012ah . and warrants deeper observations to de¬ 
termine whether gas traced by CS (s uggested to be a tr acer of 
some of the densest molecular cores: lBaan et al.ll2008l) is also 
taking part in the molecular outflow. 

To derive the dense gas column along the line-of-sight to 
the AGN in Mrk 231 (white cross on Fig. [3] moment maps), 
we And that the pixel associated with the radio point source 
contains intensity of 51.8 ± 2.OK km s“^ (using a Kelvin 
per Jansky factor of 78.35). Converting this intensity into a 
column density of H 2 (assuming an excitation temperature of 
70 K and a CS/H 2 abundance of 10“®), we And a CS-derived 
line-of-sight column to the AGN of 1/2 x Nh = A^(H 2 ) Ri 
1.1 X 10®^cm“^ (if the AGN were sitting beneath the maxi¬ 
mum value in the CS(2-1) map, the obscuring column would 
be 1V(H2) R:! 1.5 X 10®'^cm“^). While this value was consis- 
tent with a previou s estimate from X-ray absorption models 
(iBraito et al.ll2004l) . newer results from NuSTAR show Nh to 
be 1.2 X 10®^ cm“® (iTeng et al.ll2014^ . The discrepancy be¬ 
tween the CS-derived column might suggest that the obscur¬ 
ing medium in the nucleus of Mrk 231 is quite clumpy (at a 
scale much smaller than the CA RMA resolution), which is 
consistent with the suggestion by ITeng et akl (120141) that the 
wind punching out of the system has created a preferential 
line-of-sight to the AGN. 

If we were to use the CS to investigate the star formation 
efficiency within Mrk 231, the total dense gas-traced surface 
density is Eh 2 ~ 6450 Mq pc“®. Comparing this to the 
star formation surface density (assuming the dense gas and 
star formation are co-spatial) of 61 Mq yr“^ kpc“^, we And 
that this matches the expectation from the Kennicutt-Schmidt 

^ Not including t he uncertainty associated with converting CS to M(H 2 ) 
(Alatalo et alJ20r5l) 
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Fig. 1.— Integrated CO(l-O) spectrum of Mrk231 from CARMA. The 
shaded regions correspond to the ranges summed to search for wing emission, 
of 400 < V < 800km s~^ for the red-shifted wing and —800 < v < 
—400 km s~^ for the blue-shifted wing. The momentO map is shown as an 
inset in the top left, where a white cross pinpoints the location of the AGN. 

relation (lKennicuttjll998h . 

3.2. Resolved CO(l-O) outflow 

Accounting for the distance to Mrk231, 0.49" corresponds 
to a projected separation between the CO(l-O) centroids of 
415 pc. The CARMA data also successfully resolve the 
lobes, with sizes of 2.2" and 2.3" of the red- and blue- 
shifted wings, respectively. To deproject this distance, an 
accurate inclination angle is needed. Inclination angles de¬ 
rived for the molecular disk (which we assume is perpendic¬ 
ular to the ou tflow) has derived inclinat i ons angles ranging 
from i = 10° (iDow nes & S olomonl[T998HDavies et al.ir2004l) 
to i = 45° (iRichards et akl 120051) . meaning that the depro- 
jected distance ranges between 570 and 2390pc. If we as- 
sume that the outflow is conical with a large opening angle 
(iRupke & Veilleu?! 1201 U ICicone et al.l 1201^ . then the dis¬ 
tance between the wings is likely similar to the diameter of the 
lobes, which we estimate to be 1800 pc (after de-convolving 
the lobes). This size is slightly large r than the CO(l-O) sep¬ 
aration found bv ICicone et al.l(l2012h . although we see exten¬ 
sion in the East-West direction and opposite of the rotation 
of the CS disk, rathe r than the North-So uth direction as was 
observed in HCN bv lAalto et al.l (1201 2al) . While this is con¬ 
sistent with the angle of rotation of the molecular disk, the 
velocities probed by the CARMA observations exceed the es¬ 
cape velocity, and therefore are unlikely to be due to rotation. 

We use 600km s“^ for the characteristic velocity of each 
wing (as well as the de-projected radius, 900 pc) to calculate 
the dynamical timescale of the outflow of Tdyn ~ 1.5 Myr, 
about twice as long as the timescale originally reported by 
iFeruglio et al.l (120101) . If we use this d ynamical time s cale i n 
conjunction with t he outflow mass fromiFeruglio et al.l (l2010t) . 
and supported bv lAalto et al.l (1201 2al) of 5.8x10® Mq, we 
find that the total mass outflow ra t e is 3 90 Mq yr“^, a fac¬ 
tor of 2 smaller than iFeruglio et al] (1201 Oh . and a factor of ft! 2 
larger than the star formation rate in the system (IVeilleux et al.l 
l2009h . If all 390 Mq yr“^ taking part in the molecular out¬ 
flow were to escape the galaxy, we would expect the molec¬ 


ular gas to be depleted in ss 46 Myr through the action of the 
molecular outflow alone, and 32 Myr if both outflow and the 
star formation gas consumptions are considered. This deple¬ 
tion timescale assumes that all molecular gas in the thin disk is 
intercepted by the molecular outflow, which is assumed to be 
traveling perpendicular to the disk. If the molecular outflow is 
unable to interact with most of the gas in the molecular disk, 
this could increase the molecular outflow depletion timescale 
considerably. 

3.3. The importance of competing depletion mechanisms 

In Mrk231, the molecular gas is being depleted by three 
competing mechanisms: consumption through star formation, 
escape from the galactic potential via a molecular outflow and 
accretion onto the supermassive black hole. The timescale 
for star formation to completely consume the central gas is 
110 Myr. The molecular outflow appears to be the dominant 
depleter within the Mrk231 system (assuming that the out¬ 
flowing molecular mass is being completely expelled). 

While Mout of Mrk231, updated in the prior section is in¬ 
deed larger than Msfr, Mo^t might not be an accurate repre¬ 
sentation of the mass of molecular gas that is currently be¬ 
ing expelled from the system. Mrk231 is a massive sys- 
tem, with an estimated stellar mass of « 3 x 10^^ Mq 
(lU et alJ[2012h . and thus much of the outflowing molecular 
gas will not escape the potential well of Mrk231, instead 
falling back and “stirring up” the massive molecular disk in 
the center. Using a dynamical mass within 1100pc of the 
center of Mrk231 of 3.1 5 x 10^° Mq (IDownes & Solomc^ 
ll998l:lDavies et a^]l2004^ . and derive an escape velocity from 
the central kiloparsec of ~ 500km s“^. If we assume that 
all gas that is traveling above Vesc will be able to successfully 
escape the system, then 15% of th e total outflowing mo lec- 
ular gas (from the gaussian fit of IFeruglio et al.l [20101) es- 
capefl We use our updated maps to calculate the total escap¬ 
ing mass from the blue- and red-shifted wing fluxes reported 
in Section |2] and Table [H using a conserv ative merger-based 
conversion factor (iNaravanan et aljr201 Ih . and find Mgsc = 
(2.93 ± 0.49) X 10® Mq. Dividing by the dynamical time, 
we find Mesc ~ 195 Mq yr“^, which is quite comparable to 
the current star formation rate of 170 Mq yr”^ than were we 
to assume all outflowing mass is escaping. This updates the 
depletion timescale due to the escape of molecular gas, Tesc 
to 92 Myr. The combined depletion timescale rsF-i-esc is then 
49 Myr. 

The accretion onto the black hole can be derived using the 
bolometric luminos i ty of the AGN (of ~ 2.8 x 10^^ Lq\ 
IVeilleux et al.ll20(M 120131) . and assuming an efficiency rj = 
0.17, to be Macc ~ 1 Mq yr“^, or at a rate ~ 10“^ other de¬ 
pletion rates, a minor constituent to the depletion of the cen¬ 
tral molecular gas, and requiring a Hubble time to completely 
remoye the molecular gas from the system. 

Without a constant driying mechanism, the majority of the 
molecular outflow will eyentually fall back into the center, ei¬ 
ther to be accreted, formed into stars, or re-launched. One 
would expect that this gas would fall back and re-inject en¬ 
ergy into the ce ntral molecular disk, possibly acting to in hibit 
star formation (iGuillard et alJl20f5l : lAlatalo et al.ll2015h . and 

'* TFiis estimate only accounts for t[ie tlie molecular gas traced by CO(l— 
0), although the Mrk231 outflow has been detecte d in additional hace rs of 
netural gas, including OH (Sturm et alJ201 ID . NaD (Veilleux et alJ2009ll and 
H UTengetalJ2013() . 
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Fig. 2.— (Left): 6 -band archival image of Mrk231 from HST (grayscale, IKim et 3111201311 is overlaid with the integrated CARMA maps of the blue-shifted 
(—800 < V < —400 km s“^; blue) and red-shifted (400 < v < 800 km s“^; red). Contours for the red lobe is [3,4,5,6 ,7](t and [3,4,5,6,7,8 ](t for the blue 
lobe, where cr is the root mean square noise of the map. The blue and red crosses represent the centroids on the blue- and red-shifted wings, respectively. The 
white cross represents the position of the centroid of the 3mm continuum point source. Separation between the lobes is ^ 0.49", a well-resolved distance, where 
the size of the red and blue crosses represent the centroiding errors for each wing. The separation is detected to 5.8 times the centroiding error. It is of note that 
the CO(l-O) wings velocity range almost entirely exceeds the escape velocity out of the central kiloparsec of Mrk231. (Right): Integrated CO(l-O) spectrum 
of Mrk231 using emission exceeding a signal-to-noise of 3 in red- and blue-shifted wing images to create the extraction apertures. This seems to indicate that 
CARMA has successfully detected emission in the red- and blue-shifted wings. 

therefore extend the molecular gas consumption timescale. 

An extension of the depletion timescale expected for Mrk 231 
from ~ 20 Myr (the typical lifetime for O-stars) to ^ 50 Myr, 
appears more consi stent with the poststarburs t population 
found in Mrk 231 (iCanalizo & StocktonI I2000t). and com¬ 
mon i n other quasars ^Cales et al.ll201 lUCanalizo & StocktonI 
IMl, as well as other hosts of mole cular outflows (such as 
NGC 1266: lAlatalo et al.llWIll2()Il) . 

The rapidity at which Mrk 231 is currently expelling its 
molecular material compared with weaker AGNs might be 
due to the different circumstances that created this system 
and outflow. A major merger occurred in a system already 
sufficiently massive to host a massive black hole, which is 
thus energetically capable of driving the bulk of the molec¬ 
ular gas from the system. Sub-Uesc gas would be able to fall 
back into the system and inject turbulence, in effect, temporar¬ 
ily stalling efficient star formation and possibly extending the 
depletion time just slightly. This would account for the mis¬ 
match in depletion timescales between ^ 30 Myr (the deple¬ 
tion timescale should all outflowing gas escape and with the 
current star formation rat e) and ^ 50 Myr (the stellar popula¬ 
tion age seen in Mrk231: ICanalizo & Stocktonll2000h . 

Deeper observations of the molecular outflow in Mrk 231, 
focusing on spatial resolution on at different distances away 
from the AGN will be able to determine what fraction of 
the molecular gas is escaping the gravitational potential, 
and provide a more complete understanding of the driving 
mechanism as well as the energy injection rate. Low fre- 
quency imaging could pin point the location of fossil shells 
(ISchoenmakers et d]l200(lll . as Mrk 231 has already been de¬ 


tected at low frequency (ICohen et al.ll2007l) . Combining these 
observations with the multiwavelength suite available for 
Mrk 231 will provide additional constraints, such as whether 
the molecular outflow changes phase as it accelerated, and 
whether there is any evidence that gas ultimately falls back 
into the molecular disk of Mrk 231. 

4. CONCLUSIONS 

We have presented 3mm continuum, and CO(l-O) and 
CS(2-1) high resolution molecular gas maps of Mrk 231 from 
CARMA. The CS(2-1) data paint a picture of the dense 
molecular gas in Mrk 231 that is consistent with other dense 
gas tracers (lAalto et al]l2012ah . CS(2-1) predicts a dense gas 
mass (of Mdense = 1.8 ± 0.3 X 10^° Mq), and a line-of-sight 
column toward the AGN of A^(H 2 ) ~ 10^^ inc onsis- 

tent with the X-ray derived column (iTeng et alJl2014^ . pos¬ 
sibly indicative of a clumpy ISM, and the AGN-driven winds 
punching a preferential line-of-sight to the AGN. If we use the 
dense gas to derive a star formation efficiency, we And that 
Mrk 231 is consistent with the prediction of the Kennicutt- 
Schmidt relation. 

CARMA was able to resolve the individual CO(l-O) broad 
line-wings, with an inferred separation between the wings to 
be Ri 1.8 kpc given the angular size of the lobes and possi¬ 
ble inclination angles of the outflow. This is s lightly larger 
than t he CO(l-O) lobe separation reported in ICicone et alJ 
(l2012h . If we assume that not all molecular gas will suc¬ 
cessfully escape the galaxy, we can re-calculate the deple¬ 
tion timescale for Mrk231 based on the mass escape rate, of 
Ri 195 Mq yr“^. A depletion timescale using Mesc + -^sfr, 
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Fig. 3.— (Top left & middle): The integrated intensity (momentO) and mean velocity (momentl) maps of CS(2-1) within Mrk231 from CARMA. While 
there does appear to be some ordered motion, overall the spectrum appears quite turbulent. Contour levels are [0.01,0.1,0.3,0.5,0.7,0.9] times the maximum value 
of the momentO map. (Top right): The integrated CS(2-1) spectrum, with the RMS per channel shown as a small error bar on the left-hand side of the spectrum. 
The CS(2-1) line appears to have a width of Ri300 km s“^. (Bottom): The CS(2-1) channel maps. Contours are [-2.5,2.5,3.5,4.5 ]ct, where a is the rRMS of the 
CS(2-1) cube. The negative contours are gray. 


rather than Mout + -^sfr is more consistent with the post- 
starburst stellar populations that are seen in Mrk231 as well 
as many other quasars. 
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